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Abstract. The unpolarised differential cross section and the two deuteron tensor analysing powers A xx and 
A yy of the dp — > {pp}„n charge-exchange reaction have been measured with the ANKE spectrometer at the 
COSY storage ring. Using deuteron beams with energies 1.2, 1.6, 1.8, and 2.27 GeV, data were obtained 
' for small momentum transfers to a {pp} 3 system with low excitation energy. The results at the three 

lower energies are consistent with impulse approximation predictions based upon the current knowledge 
of the neutron-proton amplitudes. However, at 2.27 GeV, where these amplitudes are far more uncertain, 
agreement requires a reduction in the overall double-spin-flip contribution, with an especially significant 
effect in the longitudinal direction. These conclusions are supported by measurements of the deuteron- 
proton spin-correlation parameters C XtX and C y , y that were carried out in the dp— s- {pp} s n reaction at 1.2 
and 2.27 GeV. The values obtained for the proton analysing power Ay also suggest the need for a radical re- 
evaluation of the neutron-proton elastic scattering amplitudes at the higher energy. It is therefore clear that 
. such measurements can provide a valuable addition to the neutron-proton database in the charge-exchange 

region. 

PACS. 13.75.-n Hadron-induced low- and intermediate-energy reactions and scattering (energy < 10 GeV) 
- 25.45.De Deuteron breakup - 25.45.Kk Charge-exchange reactions 



1 Introduction 

An understanding of the nucleon-nucleon (NN) interac- 
tion is fundamental for the whole of nuclear and hadronic 
physics. The SAID database and analysis program |Tj have 
proved to be truly invaluable tools over many years for 
researchers working in this area. The general procedure 
adopted here is to take all the NN elastic scattering data 
in the literature in order to perform a phase shift analysis 
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up to a certain orbital angular momentum L max and use 
a theoretical model for higher L. When significant new 
experimental data become available, the value of L max 
can be increased. By assuming that the phase shifts vary 
smoothly with beam energy, predictions can be made for 
observables at a particular energy and it is in this way 
that the SAID program is most commonly used. 

Clearly any amplitude analysis can only be as good 
as the data used in its implementation. Though lots of 
proton-proton observables have been measured up to high 
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energies, there are significant gaps in our knowledge for 
beam energies Tjv > 2 GeV, especially at small angles. 
The situation is even more serious for neutron-proton scat- 
tering where there are major holes in the experimental 
database above about 1 GeV and data that do exist are 
not necessarily very well reproduced by the SAID pro- 
gram. For example, the only differential cross section data 
for large angle np scattering [2J, in the so-called charge- 
exchange region, seem to be consistently over-predicted in 
the SAID analysis. 

Extra information in the np sector can be found by 
performing experiments with a deuteron beam or a deu- 
terium target. The simplest of these involves measuring 
the ratio of the forward charge-exchange cross section of 
a neutron on a deuterium target to that on hydrogen, 



-Rnp(O) 



da{nd pnn) / dt 



da(np — > pn)/dt 



(1.1) 



9=0 



where t is the square of the four-momentum transfer be- 
tween the initial neutron and final proton. Due to the 
Pauli principle, when the two final neutrons are in a rela- 
tive iS-wave their spins must be antiparallel and the sys- 
tem is in the ^So state. Under such circumstances the 
nd — > p{nn} reaction involves a spin flip from the S = 1 
of the deuteron to the S = of the dineutron. In impulse 
approximation the ratio determines 



R np (0) — — 



3 \ a s f + a ns f 



(1.2) 



where a s f (ens/) is the forward np charge-exchange differ- 
ential cross section with (without) a spin- flip [3J. 

Extensive measurements of R np (0) have been reported 
and these have recently been extended up to 2 GeV at 
Dubna pQ. Although impulse approximation predictions 
for R np (0) on the basis of the current SAID amplitudes re- 
produce well these observations from 300 MeV up to about 
800 MeV, there are serious discrepancies above 1 GeV [5] . 
In terms of Eq. (|1.2p . it seems that the SAID solution 
overestimates the contribution of the spin-flip amplitudes 
to the forward np charge- exchange cross section. 

More detailed information on the np charge- exchange 
amplitudes can be derived by using a polarised deuteron 
beam or target and studying the dp — > {pp} s n reaction [BJ . 
To achieve maximum sensitivity, the excitation energy E pp 
in the final pp system must be very low so that the final 
diproton is in the ^Sq configuration. Experiments from a 
few hundred MeV up to 2 GeV [3[8] have generally borne 
out well predictions based upon the impulse approxima- 
tion [BJ. 

In order to constrain the np amplitudes using this 
approach, the ANKE collaboration has embarked on a 
systematic programme to measure the dp — > {pp} s n ob- 
servables up to the maximum COSY deuteron energy of 
Td s» 2.3 GeV [3]. The proof of principle of the method 
was the experiment carried out at a deuteron energy of 
1.17 GeV where, because of the wealth of neutron-proton 
data, the SAID amplitudes used as input in the calcula- 
tions should be quite reliable. The measured values of the 



unpolarised cross section and the two deuteron Cartesian 
tensor analysing powers A xx and A yy were then quantita- 
tively reproduced in impulse approximation [10] . Dilutions 
of the signals due to higher partial waves in the final pp 
system were taken into account in the calculations jll) . 

The results of similar measurements are presented here 
at deuteron beam energies of 1.2, 1.6, 1.8, and 2.27 GeV. 
Whereas the unpolarised differential cross section is cor- 
rectly predicted in impulse approximation at the three 
lower energies, this gives a result that is about 15% too 
high at 2.27 GeV. Such a factor is in fact consistent with 
the discrepancy between the SAID np predictions and the 
data in the charge-exchange region [2] , taking into account 
the SAID overestimation of the spin-flip contributions, as 
measured in the inclusive R np (0) experiment [H[S]. 

The situation is rather similar for A xx and A yy , where 
impulse approximation reproduces the measurements very 
well at 1.8 GeV and below but gives a much poorer de- 
scription at 2.27 GeV. Agreement can be restored if the 
spin-longitudinal amplitude is reduced by ~ 25% com- 
pared to the SAID values. This is not entirely unexpected 
because the deuteron tensor analysing power in the for- 
ward direction is closely linked to the longitudinal and 
transverse spin-transfer coefficients in np charge exchange. 
Since there are no measurements of these parameters at 
small angles in the T n ps 1.135 GeV region, the SAID am- 
plitudes are fixed here mainly by the assumptions in the 
model. 

Some confirmation of the renormalisation of the spin- 
longitudinal amplitude is offered by the values of the two 
deuteron-proton spin-correlation parameters C x , x and C VyV 
that were measured at 1.2 and 2.27 GeV by replacing the 
hydrogen cluster target with a specially designed cell that 
can be filled with polarised gas [12] . In the So limit these 
parameters depend upon interferences between the spin- 
spin amplitudes and, although the SAID input reproduces 
well the 1.2 GeV data, at 2.27 GeV a 25% reduction in the 
spin-longitudinal term leads to a much better description. 

The polarised hydrogen target [13| also allowed us to 
measure the proton analysing power A v y in the dp '— > {pp} s n 
reaction and the conclusions to be drawn here are broadly 
similar. Despite the AP y signal being quite small over the 
ANKE angular range, it is reproduced quantitatively in 
the impulse approximation calculations at 1.2 GeV, though 
these fail badly at 2.27 GeV. This observable depends 
upon the interference of a spin-orbit amplitude with a 
spin-spin term and it would seem that the current SAID 
analysis [I] grossly underestimates the imaginary part of 
the spin-orbit amplitude at the higher energy. 

We have described the phenomenology of the deuteron 
charge-exchange reaction at some length in earlier publi- 
cations 10,14,15 but, in order to make the current pa- 
per more self-contained, some of this material is repeated 
in sect. [2] The experimental set-up for measuring with 
the hydrogen cluster-jet target is identical to that used 
in our earlier work |10) and so most of the emphasis in 
sect. [3] is on the polarised hydrogen gas cell used in the 
spin-correlation measurements. The measurements of the 
luminosity through the observation of the quasi-free pn — > 
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dir° reaction are reported in sect. |TJ where the results for 
the unpolarised dp — > {pp} s n differential cross sections are 
compared with impulse approximation calculations. 

The polarisation of the deuteron beam was established 
at 1.2 GeV and, since there are no depolarising resonances 
for the deuteron in the COSY energy range, the analysing 
powers for the dp — > {pp} s n reaction could be measured at 
various energies and the results are presented in sect. [SJ 
The polarisation measurements with the gas cell have no 
parallel in our earlier work. The method used here relies 
on the data taken with the polarised deuterium target 
and the measurement of the analysing power of quasi-free 
np — > dir°. Taken together with the measurement of the 
vector polarisation of the deuteron beam, this allowed us 
to extract the spin-correlation and proton analysing power 
results given in sect. [5] Our conclusions regarding the use- 
fulness of charge exchange on the deuteron in the study 
of neutron-proton elastic scattering amplitudes are drawn 
in sect. [7] 



2 Deuteron charge exchange in impulse 
approximation 

The cross section and spin observables for the dp — > {pp} s n 
reaction have been extensively discussed in the literature [51 
ITT] and only essential formulae are collected here. In im- 
pulse approximation the deuteron charge exchange ampli- 
tude is proportional to a neutron-proton charge exchange 
amplitude times a form factor that reflects the overlap be- 
tween the initial deuteron wave function and that of the 
outgoing diproton system. 

The elementary np — > pn amplitude may be written in 
terms of five scalar amplitudes in the cm system as 

fn P = a(q) + i^(q)((T 1 + <x 2 ) • n + ,3(g)(cr 1 ■ n)(er 2 • n) 
+6{q){<r 1 • m)(<r 2 • m) + e{q){a 1 ■ l){a 2 ■ I), (2.1) 



where q = \f—i is the three-momentum transfer and the 
Pauli matrices <x are sandwiched between neutron and 
proton spinors. Here a is the spin-independent amplitude 
between the initial neutron and final proton, 7 is a spin- 
orbit contribution, and /?, S, and e are three spin-spin 
terms. 

The orthogonal unit vectors used in Eq. (I2.1[) are de- 
fined in terms of the initial neutron (K) and final proton 
(K') cm momenta; 



K x K' 

n ~ \K x K'\ ' 



m = 



K' K 

\K'-K\ 



I = 



K' + K 
\K' + K\ 



(2.2) 



The amplitudes are normalised such that the 
differential cross section has the form 



np 



pn 



da 

df 



\a(q)\ 2 + \P(q)\ 2 +2h(q)\ 2 + \S(<l)\ 2 + m\' 



np—tpn 



(2.3) 

In the limit of very low E pp = k 2 /m p , where k is 
the pp relative momentum, the deuteron charge exchange 



necessarily involves a spin flip from the S = 1 of the 
deuteron to the S = of the diproton. In this case the 
contribution from the spin-independent amplitude a drops 
out and one is left with only the spin-flip cross sections, 
as in the sum rule of Eq. (I1.2[) . The observables involv- 
ing only the initial spins that are accessible at ANKE are 
linked to the amplitudes through [TT1I16] : 



dkPk 
I A* 



±l{S~(k,±q)} 2 , 



= 



IA% = -27m(/3* 7 ) 



T A — 

T A — 

I C x ,x = 

Cyy,y = 



l/3| 2 + M 2 + H 2 
\S\ 2 R 2 + \e\ 2 -2\ 
-2Re(e*5)R, 
-2Re(e*P), 



■2\6\ 2 R 2 , 
| 2 -2| 7 | 2 



where the spin-flip intensity 
I = 



2 + l7| 2 + H 2 



\6\ 2 R 2 . 



The function 



R= S+(k,±q)/S-(k,±q) 



(2.4) 



(2.5) 



(2.6) 



is the ratio of two transition form factors that involve the 
S- and Z?-states of the deuteron wave function. In the 
forward direction R = 1. 

Since 7(g) vanishes in the forward direction, the con- 
tributions of |7((?) | 2 to the cross section and the Cartesian 
tensor analysing powers A xx and A yy in Eq. (|2.4[) are al- 
most negligible under the conditions of the ANKE exper- 
iment. Measurements of the unpolarised cross section and 
the two transverse analysing powers can therefore deter- 
mine separately the values of |/3(g)| 2 , |e(<7)| 2 , and |<5(q)| 2 at 
fixed momentum transfer q. The two deuteron-proton spin 
correlations that are measurable at ANKE, C X:X and C VtVl 
fix two of the relative phases. The proton analysing power 
A v y gives mainly information on the spin-orbit amplitude 

7(9)- 

Although the formulae given here describe the general 
features of our data, detailed comparisons with theory are 
made using a program that takes higher final pp waves into 
account. These can, in particular, dilute the polarisation 
signals [TT] . 



3 The experimental facility 



The experiments reported here were carried out over three 
different time periods using the ANKE magnetic spec- 
trometer [17] that is placed at an internal target position 
of the COoler SYnchrotron (COSY) [T5] of the Forschungs- 
zentrum Jiilich. Initially a polarised deuteron beam was 
used in conjunction with an unpolarised hydrogen cluster 
target [JJ5]. In 2005 the dp — > {pp] s n reaction was studied 
at deuteron beam energies Td — 1.2, 1.6 and 1.8 GeV. 
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The following year, the beam energy was increased to 
2.27 GeV, with 1.2 GeV being repeated for polarimetry 
purposes. The equipment used was described in our earlier 
publications PUIH^IITS] . However, for the study of the spin- 
correlation parameters in 2009 [2U] , a newly developed po- 
larised internal target [T2"UT5] was installed at ANKE and 
this was employed in experiments at 1.2 and 2.27 GeV. 

3.1 The polarised deuteron beam at COSY 

The polarised deuterium ion source at COSY provides 
beams with different spin configuration (15) . It uses ra- 
dio frequency transition units and quadrupole magnets to 
exchange the occupation numbers of the hyperfine states 
in the atom. The source was set up to provide a variety of 
states with different tensor and vector polarisations but, 
as listed in Table [U the selections differed for the experi- 
ments with cluster or polarised cell target. 

The COSY cycles were configured to provide beam first 
at 1.2 GeV and then, without additional injection, acceler- 
ate the deuterons to one of the higher energies. This proce- 
dure allows the use of the polarisation export method [2T] , 
which is crucial in the measurement of spin observables at 
higher energies. This technique involves undertaking the 
polarimetry measurements at the lowest Td = 1.2 GeV 
flat top energy, where the analysing powers are precisely 
known, and assuming that the beam polarisation is un- 
changed at the higher energy. This procedure is viable be- 
cause there are no depolarising resonances for deuterons 
in the COSY energy range. This was checked at the 4% 
level by repeating the measurement of the analysing pow- 
ers after deceleration [2Tj . 

3.2 The polarised hydrogen cell target at ANKE 

The ANKE polarised internal gas target [T31I22] uses an 
Atomic Beam Source (ABS) that is capable of producing 
both polarised hydrogen and deuterium beams [23] - The 
first tests with the hydrogen atoms gave polarisations of 
+0.89 ±0.01 and —0.96 ±0.01 for spin-up and spin-down, 
respectively [T51l2~4"] . 

The polarised atomic beam could be used directly from 
the source as a jet target. During the commissioning runs, 
the ABS demonstrated an integral jet-target thickness of 
about 1.5 x 10 11 cm~ 2 , which is consistent with the pre- 
dicted value 24 . However, much higher target densities 
can be achieved if one uses a storage cell fed by the ABS. In 
order to achieve the maximum density, it is important to 
minimise the dimensions of the storage-cell tube. But, on 
the other hand, this limits the number of particles stored 
because of the beam heating and the consequent losses 
on the cell walls. A maximum target density of about 
10 13 cm~ 2 was achieved with the cell during the com- 
missioning runs and this resulted in luminosities of up to 
10 29 cm _2 s _1 , depending upon the beam intensity [22"] . 

During the 2009 beam time a cell made of 25 /j,m 
thick aluminium foil (99.95% Al) was used [12] . In order 
to minimise depolarisation on the cell surface, its inner 
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Table 1. The different configurations of the polarised deuteron 
ion source used in experiments carried out with the unpolarised 
cluster target, showing the nominal (ideal) values of the vec- 
tor (P z ) and tensor (P zz ) polarisations and relative beam in- 
tensities. For the spin-correlation measurements, the more re- 
stricted set was used. Note, that P z and P zz are labeled con- 
ventionally in the reference frame of the source [15] . 



walls were coated with Teflon. The cell had dimensions 
VxVxZ = 20xl5 x 370 mm 3 , where Z is measured 
along the beam direction with X and Y referring to the 
horizontal and vertical transverse directions, respectively. 

A dedicated beam development was required to ensure 
that the COSY beam passed successfully through the cell. 
Electron cooling [25] and stacking injection [27] were em- 
ployed, with hundreds of injections per cycle to increase 
the number of stored deuterons in the beam. In order to 
avoid excessive background coming from the interactions 
of the beam halo particles with the cell wall, scrapers were 
installed upstream of the target region. 

In the hydrogen case, the ABS was configured to pro- 
duce two polarised states with equal gas densities. The 
polarisation of the hydrogen target was flipped between 
spin-up (t) and spin-down (!) every five seconds through- 
out the whole COSY cycle, which lasted for one hour. 
Such a procedure simplifies the later analyses by obviat- 
ing the need to consider the luminosities while calculating 
the asymmetries between these states. In order to have 
a possibility to measure the target spin-up f and spin- 
down I polarisations separately, runs with an unpolarised 
hydrogen cell-target were also undertaken. 

3.3 The ANKE detection system 

The forward detector system (FD) of the ANKE magnetic 
spectrometer, illustrated in Fig.Q] is used for the deuteron 
charge-exchange studies. The multiwire chambers in the 
FD serve for track reconstruction and the three layers of 
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MWDC MWPCs 



s Hodoscope 



Fig. 1. The ANKE experimental set-up showing the positions 
of the three dipole magnets Dl, D2, and D3. The Forward De- 
tector (FD) consists of multiwire drift (MWDC) and multiwire 
proportional (MWPC) chambers, and a hodoscope composed 
of three layers of scintillation counters. For measurements with 
the polarised target, the ABS replaced the cluster-jet in the 
target chamber, with the detection system being identical in 
the two cases. The axes of the coordinate system are indicated. 



the scintillation hodoscope permit the measurements of 
the arrival time and energy-loss that are required for par- 
ticle identification |28j . The main trigger used in the ex- 
periments consisted of a coincidence between the different 
layers in the hodoscope of the FD. Figure [2] shows the ex- 
perimental yield of ANKE for single charged particles at 
Td = 1.2 GeV in terms of the laboratory production an- 
gle in the horizontal plane and the magnetic rigidity. The 
kinematic curves for some of the possible nuclear reactions 
are also illustrated. 
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Fig. 2. Scatter plot of singly charged particles detected in 
ANKE from the interaction of 1.2 GeV deuterons with a hy- 
drogen cluster-jet target in terms of the laboratory production 
angle in the horizontal plane and the magnetic rigidity. The 
loci corresponding to four common nuclear reactions are also 
shown. These include the dp —¥ {pp} a n reaction at zero E pp . 



dp — > {pp} s n and the quasi-free dp — > p sp dir°, where the 
proton, p sp , has about half the beam momentum. The 
latter reaction is used to measure both the vector po- 
larisation of the deuteron beam or hydrogen target and 
also the luminosity. After recording two charged parti- 
cles, deuteron-proton pairs are separated from the remain- 
ing two-track events (mainly proton pairs) in the subse- 
quent analysis by using the time information from the ho- 
doscope. As demonstrated in Fig. [31 if one assumes that 
both detected particles in the pair are protons, the calcu- 
lated (AT C ) time of flight difference from the target does 
not match with the measured one (AT m ) for other pairs. 
After recognising the two charged particles, the missing- 
mass distribution allows one to identify the reaction. 
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Fig. 3. Scatter plot of measured (AT m ) and calculated (AT C ) 
time differences between pairs of charged particle registered in 
the ANKE forward detector at T d = 1.2 GeV. AT C was calcu- 
lated assuming that both particles were protons. The separa- 
tion between dp and pp pairs is very clear. 



4 The cross section determination 

The cross section a for a given physical process is given 
in terms of the corresponding counting rate R and the 
luminosity L through: 



a = R/L 



(4.1) 



The luminosity, which is the product of the target den- 
sity and beam intensity, can be measured in various ways. 
In the current analysis we relied on the measurement in 
parallel of the rate for a process with a well-known and 
sizeable cross section. Once the luminosity is known, ab- 
solute values of cross sections for other reactions can be 
deduced from the count rates measured in the experiment. 



4.1 The dp — > p sp dTv° reaction 



Among the reactions observed, there are two that are The dp — » p sp dir° reaction is used in the determination of 
of particular interest, namely the deuteron charge-exchange the luminosity in this experiment. This is identified in the 
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ANKE forward detector by detecting both charged par- 
ticles (cf. Figs. [2] and [3]). After recognising the dp-pairs, 
the reaction is finally isolated on the basis of the missing- 
mass distributions [10 . There is an accidental background 
at very small \AT C \ that is randomly distributed in AT m . 
This is caused by fast particles, mainly protons, that are 
produced in a different beam-target interaction. The con- 
tribution from such accidental events in the vicinity of the 
fast dcuteron branch of the dp — > p sv dir° reaction increases 
rapidly with energy. It varies between 18% and 30% for 
T d > 1.6 GeV, whereas it is less then 3% at 1.2 GeV. 
The background is negligible at all energies for the slow 
deuteron branch. 

The properties of background were studied using the 
data for which \AT C \ < 2 ns and \AT m \ > 12 ns. As 
can be seen from Fig. [3] no true coincidence two-track 
events are expected in this region. These data provided 
the shape of background in the distributions of the miss- 
ing mass and the deuteron laboratory scattering angle in 
the dp — > p S pdir° reaction. The normalisation for the back- 
ground was found by comparing the background missing- 
mass spectra with that for the identified dp pairs. The nor- 
malised background was then subtracted from the dp — > 
p S pdir° angular distributions. 




P p [GeV/c] 

Fig. 4. The momentum distribution of the fast proton from 
the dp — >• p Bp dir° reaction at Td — 1.2 GeV, transformed into 
the rest frame of the incident deuteron is compared with the 
Monte Carlo simulation (solid histogram). 



At intermediate energies, soft deuteron collisions are 
generally dominated by the interaction of one of the nu- 
cleons in the nucleus, the other nucleon being a spectator. 
When the proton acts as a spectator, p sp , the dp —} p sp dTr° 
reaction can be interpreted in terms of quasi-free np — > 
dir° pion production. To confirm the spectator hypothe- 
sis, a Monte Carlo simulation has been performed within 
PLUTO [inj using the Fermi momentum distribution from 
the Paris deuteron wave function |30j . As is clear from 
Fig. HI the data are consistent with quasi-free production 
on the neutron leading to a spectator proton. However, in 
order to reduce further possible contributions from mul- 
tiple scattering and other mechanisms, only events below 
60 MeV/c were retained for the luminosity evaluation. 



The determination of the angles for the quasi-free np — > 
dw° reaction is complicated by the Fermi motion of the 
nucleons inside the deuteron. Due to this effect, the effec- 
tive neutron beam energy, T„, is spread around half the 
deuteron beam energy with a width arising from the Fermi 
momentum. At a beam energy of 600 MeV per nucleon, 
the FWHM is 90 MeV for a p sp < 60 MeV/c cut. Further- 
more, the neutron direction is not precisely aligned along 
that of the beam, but is spread over some solid angle. Since 
this introduces an incident angle, which is several degrees 
in the laboratory system (depending on the beam energy), 
it has to be taken into account. These considerations ap- 
ply to both the polar and azimuthal angles. In order to 
correct for this effect, the three- momentum of the incident 
neutron was reconstructed using the information from the 
spectator-proton momentum. The deuteron polar angle 
was measured from the neutron momentum instead of the 
beam direction. The azimuthal angle was defined between 
the normals to the COSY ring and deuteron scattering 
plane. 

Isospin invariance requires the cross section for np — > 
dir to be half of that for pp — > dir + , for which there are 
numerous measurements |31j . An additional advantage of 
using this reaction for normalisation is that the typical 
5% shadowing effect in the deuteron (where one nucleon 
hides behind the other) should be broadly similar in the 
dp — > {pp}X and dp — > p sp dTr a reactions. 

In order to investigate the acceptance of the ANKE 
forward detector for different reactions, a full simulation 
was performed based on GEANT software [32] • The same 
track reconstruction algorithm was used in the simula- 
tion and the data analysis. In order to get as precise a 
description of the experiment as possible, the dispersion 
of the hits in the MWPC, the background hits produced 
by accidental coincidences, and the noise in the multiwire 
chambers readout electronics, as obtained from the experi- 
mental data, were also included in the simulation [25] ■ The 
quality of the simulation may be judged from the distribu- 
tions of the deuteron production angle in the laboratory 
system that are shown in Fig. [S] 

4.2 Luminosity measurements 

To determine the luminosity we use Eq. (I4.1[) and insert 
detailed expressions for R and a: 

R -^cxp (A/tot/Wacc) lA s 

where R cxp is the count rate from the quasi-free np — > 
dir° reaction, corrected for the trigger dead time. A acc 
is the number of counts in the simulation that pass all 
the criteria used in the experimental data processing and 
A^tot is the total number of simulated events. These are 
summed over the neutron kinetic energy, subject to the 
p sp < 60 MeV/c cut, and over the given angular range. 
The da($, T n )/ dd differential cross section was taken from 
the pp — > d-K + database [31] , where $ is the deuteron polar 
angle in the laboratory frame. 
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4.3 The dp — > {pp} s n cross section 

Having identified two fast protons in the final state and 
selected low E pp events, the dp — > {pp} s n reaction was iso- 
lated on the basis of the missing-mass distributions, which 
are shown at three energies in Fig. |51 In addition to the 
dominant neutron peak, there are also many events for 
M x > 1080 MeV/c 2 that must correspond to pion produc- 
tion [33 . However, very few of these leak into the neutron 
region and the background from this under the neutron 
peak is at most at the per cent level. This is also true for 
data in the individual momentum-transfer q bins. Random 
background, which was studied using the timing informa- 
tion, was at the 1-3% level and could be easily subtracted. 



Fig. 5. Simulated (solid histogram) and experimental angular 
distributions for the dp —¥ p B pdiv reaction at Td = 1.2 GeV. 
The left and right panels correspond to fast and slow deuterons, 
respectively. 

The FD detector acceptance changes rapidly with an- 
gle. In order to minimise systematic errors, the total an- 
gular range was binned and the luminosity evaluated sep- 
arately for each bin. Data at the acceptance edges (the 
smallest and the largest angles) were less reliable, due to 
the greater uncertainty in the evaluation of the accep- 
tance, and showed systematic shifts in luminosity. Such 
angular intervals were discarded and the average recom- 
puted. 

The values of the average luminosities determined from 
the np — > dir° reaction are given in Table [5] The er- 
rors quoted include statistical ones from the experimen- 
tal counts and those introduced by the background sub- 
traction procedure. Uncertainties coming from the SAID 
database [31], which were estimated by studying the ex- 
perimental results in the relevant regions, are listed sepa- 
rately. 



[GeV] 


Average luminosity 
[cm-V 1 ] 


Measurement 
uncertainty 

[%] 


SAID 

uncertainty 
[%] 


1.2 


1.76 x 10 30 


1.1 


2.2 


1.6 


1.84 x 10 31 


2.0 


5.1 


1.8 


1.61 x 10 31 


2.8 


4.4 


2.27 


1.18 x 10 30 


5.0 


3.8 



Table 2. Average luminosities achieved with the cluster-jet 
target at four different beam energies. Shown separately are 
the uncertainties associated with the measurement and with 
the experimental data used as input in the estimations. 
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Fig. 6. The missing-mass M x distributions for the dp — > 
{pp} s X reaction at three deuteron beam energies. The back- 
ground under the neutron peak is a negligible. 



The dp — > {pp} s n cross section determination was per- 
formed in a similar manner to the luminosity evaluation 
that used the dp — > p sp dir° reaction. It involved the same 
technique for correcting the experimental count rates and 
estimating the detector acceptance. Since the dp — > {pp} s n 
reaction has a three-body final state, in principle the cross 
section is a function of five independent variables. How- 
ever, within the impulse approximation, by far the most 
important of these are the excitation energy E pp in the 
final pp diproton and the momentum transfer q from the 
proton to the neutron. In accordance with Eq. (I4.ip . the 
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two-dimensional differential cross section was evaluated in 
terms of the integrated luminosity Lj nt from: 



d 2 a(q,E p 
dq dEpp 



1 N eX p(q,E pp 

Lint 



N tot (q,E P p) 



N acc (q,E P p)Aq AE p 



(4.3) 



where N cxp (q, E pp ) is the corrected number of experimen- 
tal events for given values of q, measured in laboratory 
system, and E pp . N to t(q, E pp ) and N acc (q, E pp ) are the to- 
tal and accepted numbers of simulated events respectively. 
Aq and AE pp correspond to bin widths in momentum 
transfer and excitation energy, respectively. The resolu- 
tions in q and E pp are about 4 — 8 MeV/c and better than 
0.3 MeV, respectively. 

The cross sections were further integrated over E pp < 
3 MeV in order to provide the da/dq differential distribu- 
tion presented in Fig. [3 This includes also the new results 
obtained at Td = 1.2 GeV. In addition to the statisti- 
cal errors arising from the experimental count rates that 
are shown, there are also overall systematic uncertainties 
arising from the luminosity determinations, given in the 
Table [2] Within these uncertainties, the agreement with 
the theoretical impulse approximation predictions |llj at 
Td = 1.2, 1.6, and 1.8 GeV is very encouraging and is in 
line with similar data analysed at 1.17 GeV [TU]. In con- 
trast, the unpolarised differential cross section at Td = 
2.27 GeV falls about 15% below the predictions based 
upon the current np — > np partial wave analysis pQ. As 
we shall see later, similar discrepancies are found in the 
spin observables of the dp — > {pp} s n reaction but only at 
this highest energy. 



5 Determination of deuteron analysing 
powers 




150 200 
q [MeV/c] 

Fig. 7. Differential cross sections for the dp —¥ {pp} 3 n reac- 
tion at four different energies compared with impulse approx- 
imation predictions based upon the current SAID np — > np 
amplitude analysis ,1] . The data are integrated over the E vv < 
3 MeV interval. Only statistical errors are shown. 



5.1 Tensor polarisation of the deuteron beam 



In order to measure the deuteron analysing powers for the 
dp — > {pp} s n or other reaction, the first step has to be the 
identification of the polarisations of the various deuteron 
beams used in the experiment. With complete efficiencies 
in the transition units, the polarisations should approach 
the ideal values given in Table [TJ However, this is never 
the case in practice and the beam polarisation must be 
determined separately for each of the beams used. In our 
earlier work at 1.17 GeV [TU], a variety of nuclear reactions 
with known analysing powers were measured and these 
were used to establish values for the polarisations. These 
showed that the analysing powers for the dp — > {pp} s n re- 
action were well reproduced in the impulse approximation 
calculations [UJ. Since the deuteron charge-exchange re- 
action can be so well identified and measured at ANKE, 
we use this reaction itself to measure the beam tensor 
polarisation at the neighbouring energy of 1.2 GeV, the 
necessary analysing powers being taken from the impulse 
approximation estimates. Apart from the large counting 
rates, this approach has the advantage of being insensitive 
to the deuteron vector polarisation. 



The tensor polarisation of the deuteron beam was mea- 
sured using the dp — » {pp} s n reaction at Td = 1.2 GeV. 
The circulating deuteron beam was polarised perpendic- 
ularly to the horizontal plane of the machine. As already 
mentioned, the beam vector (P z ) and tensor (P 2Z ) polari- 
sations are labeled conventionally in the reference frame of 
the source. In contrast, all the spin observables discussed 
later refer to the right-handed coordinate system of the 
reaction frame, where the beam defines the Z-dircction 
while the stable spin axis of the beam points along the 
V-direction, which is perpendicular to the COSY orbit. 

No dependence is expected on the vector polarisation 
for small E pp [6] and this was checked in our earlier ex- 
periment [10 . The numbers N(q, 4>) of diprotons produced 
at momentum transfer q and azimuthal angle </> with re- 
spect to the V-direction are given in terms of the beam 
polarisation by 

= C n {1 + \P ZZ [A xx (q)(l - COS 20) + 

N {q) 

A yy (q)(l + cos2 ( f,)}}, (5.1) 
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where Nq (q) are the numbers for an unpolarised beam and 
C n is the relative luminosity of the polarised beam. 

During the course of the polarised measurements, var- 
ious configurations of the ion source were used and the 
beam polarisation had to be determined separately for 
each state. Several methods to fix the relative luminosities 
C n of a state with respect to the unpolarised mode are pos- 
sible at ANKE [21]. The one provided by the dp -> p sp X 
reaction is preferable because the number of single-track 
events is enormous for all beam energies. No dependence of 
the rates on the tensor polarisation of the deuteron beam 
was found for proton spectator momenta below 60 McV/c. 

The break-up data were divided into several bins of 
momentum transfer q and distributions in cos 20 were con- 
structed for each bin and polarisation mode. The ratios to 
the unpolarised state were fitted using Eq. (15.11) . the the- 
oretical predictions for A xx and A yy being taken at mean 
values of q in each bin. The validity of this approach was 
checked at Td = 1.17 GeV in the earlier experimental stud- 
ies at ANKE [JJ] • The beam polarisation in each state was 
taken as the weighted average over the different values of 
the momentum transfer. 

The maximum values of P zz were ~ 85% of the ideal 
values for the 2005 data. But, in the 2006 data, the max- 
imum tensor polarisation dropped to rj 55% of the ideal, 
with little change in the vector polarisation. The differ- 
ence has been ascribed to the changed efficiencies of the 
units in the COSY deuterium ion source [21] . For the high 
\P ZZ \ modes in Table [TJ the estimated polarisation var- 
ied in (0.75 — 0.85) range in the 2005 data and in the 
(0.39 — 0.60) range in the 2006 data. The typical uncer- 
tainty in P zz was 0.02 - 0.04. 

The vector analysing power of the dp — > {pp} s n reac- 
tion is predicted to vanish in the 1 Sq limit [B] and this 
was confirmed at 1.17 GeV in our earlier work [ID] . As a 
consequence, the vector polarisation of the beam is unim- 
portant for the tensor analysing power studies carried out 
with the cluster target. This is no longer the case for the 
spin-correlation measurements with the polarised cell tar- 
get. The determination of the vector polarisation of the 
deuteron beam in this case is described in sec. H) 



5.2 Deuteron tensor analysing powers 

The deuteron Cartesian tensor analysing powers A xx and 
Ayy were extracted using Eq. (|5.1[) in much the same way 
as for the polarimetry, with the beam tensor polarisation 
P zz being determined at Td — 1.2 GeV. The ratios of the 
polarised to unpolarised corrected count rates were fitted 
in terms of the two free parameters A xx and A yy . The 
procedure was repeated for different polarisation states 
and results averaged over these source modes. 

The results for the tensor analysing powers are shown 
in Fig.[S]at three beam energies as functions of the momen- 
tum transfer. The agreement between the experimental 
data and the impulse approximation predictions is very 
good at T n = 800 and 900 MeV. At these energies the 
SAID np amplitudes pQ used as input in the calculations 
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Fig. 8. Tensor analysing powers A xx (squares) and A yy (tri- 
angles) of the dp — >• {pp} s n reaction at three beam energies for 
low diproton excitation energy, E pp < 3 MeV, compared to im- 
pulse approximation predictions based upon the current SAID 
np — » np amplitudes pQ. The dashed curves at 2.27 GeV cor- 
respond to a uniform reduction of the spin-longitudinal s(q) 
amplitude by 25%. The error bars include the uncertainties 
from the beam polarisation and relative luminosity C n - 



are considered to be reliable but at 2.27 GeV, the maxi- 
mum deuteron energy available at COSY, the agreement 
is much worse. Since there are also severe discrepancies in 
the unpolarised cross section at this energy, it is natural to 
question whether there might be deficiencies in the SAID 
np analysis at this energy. 

The experimental value of A xx = A yy in the forward 
direction (q w 0) is significantly more negative than the 
predictions using the SAID np amplitudes at 1.135 GeV. 
However, it can be seen from Eq. (|2.4[) that a relative re- 
duction in the s amplitude might improve the predictions. 
To check this possibility, the predictions were recomputed 
with the e(q) amplitude being reduced uniformly by 25%. 
This gives the much better overall agreement with the 
data that is demonstrated by the dashed curves in the 
lower panel of Fig. [S] This therefore suggests that the cur- 
rent SAID amplitudes [1 might overestimate the relative 
strength of the e(q) at small q but further proof is re- 
quired and this is furnished by the measurements of the 
spin correlations. 
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6 Determination of the deuteron-proton 
spin-correlation parameters 

In experiments with the unpolarised cluster-jet target, the 
dp — > {pp} s n reaction is only sensitive to the tensor polar- 
isation of the beam and the values of this could be es- 
tablished by using the dp — > {pp} s n reaction itself at the 
1.2 GeV calibration energy, as described in sec. [SJ In con- 
trast, in order to determine the spin-correlation parame- 
ters C XtX and C ViV , one has first to determine the vector 
polarisation of the deuteron beam as well as the polarisa- 
tion of the hydrogen in the target cell. The basis of both 
measurements is the quasi-free np — > dn° reaction which, 
at small spectator momenta, is insensitive to the deuteron 
tensor polarisation. 



for each event. For a two-track event, this can be done 
with the use of the arrival-time difference for two particles, 
measured in the scintillation hodoscope. In our kinemati- 
cal conditions, where the deuteron is at least twice as slow 
as the proton, such a difference is a sensitive function of 
Z. The three-momenta of the two particles and the Y and 
Z coordinates of the vertex are found through an overall 
fit procedure that uses the information from both the wire 
chambers and the hodoscope. 

Figure [9] shows the distribution of interaction points in 
the Y-Z plane. In addition to helping in the angular de- 
termination, the vertex reconstruction allows one to make 
cuts along the cell axis to minimise the background from 
the rest gas that is spread throughout the target chamber. 



6.1 The beam and target polarimetry 

The polarisation of the hydrogen target and the vector 
polarisation of the deuteron beam were both determined 
using the quasi-free np — > dir° reaction. Now it is well 
known that, if one integrates over all Fermi momenta in- 
side the deuteron, the nucleon polarisation in the deuteron 
P™ would be reduced from that of the deuteron Pf by a 
factor 

P? = (l-lPn)P?, (6.1) 

where Pd is the deuteron D-state probability. However, 
since the D-state effects vanish like (p sp ) 2 , the dilution of 
the polarisation signal by the deuteron D-state is negli- 
gible if only data with p sp < 60 MeV/c are used in the 
subsequent analyses. Such a cut preserves a large part of 
the statistics. 

For an unpolarised deuteron beam incident on a po- 
larised hydrogen target with spin-up (f) and spin-down 
(4-), the asymmetry ratio e between polarised N 1 *^ and 
unpolarised N° yields has the form: 



iV*(+)(fl,0) 
N°(9) 



= I + Q ni) A y (9) cos (/>, (6.2) 



where 9 and <fi angles are polar and azimuthal angles, re- 
spectively, of the detected particle and Q is the target po- 
larisation. Isospin invariance requires the analysing powers 
in the np —¥ dir° and pp — > dir + reactions to be identi- 
cal and there are numerous measurements of the proton 
analysing power A y {9) of the latter in the 600 MeV re- 
gion [31] . Data in the centre-of-mass system were used in 
the analysis. The polar d and azimuthal <fi angles of the 
deuteron were determined according to the procedure de- 
scribed in sec. 14.11 



6.2 The dp — > p sp dn reaction with the cell target 

The cell introduces additional complications in the deter- 
mination of the angles because of the spread of the in- 
teraction points along the cell axis. The reconstruction of 
the longitudinal vertex coordinate Z is therefore required 
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Fig. 9. Vertex reconstruction in the YZ plane using corre- 
lated deuteron-proton pairs. The rectangle shows the physical 
dimensions of the cell (Y x Z = 15 x 370 mm 2 ). 



A second major complication arises from the scatter- 
ing of the beam halo particles on the cell walls. This 
can produce additional background that would dilute the 
analysing power signal. As mentioned earlier, the dedi- 
cated beam development enabled the bulk of the beam 
to pass through the cell without hitting the walls. For 
this reason, recording data with an empty cell would take 
much more time to collect sufficient statistics to determine 
the background. Additional runs were therefore recorded 
where nitrogen gas was injected into the cell to simulate 
the shape of the background (see details in Ref. [H]). The 
background subtraction was performed for each polarisa- 
tion state by using the missing-mass distributions for the 
hydrogen and nitrogen data. One such example is shown 
in Fig. QUI 



6.3 The target and beam polarisation 

The target polarisation was measured using data taken 
with an unpolarised deuteron beam at Td = 1.2 GeV. Af- 
ter vertex reconstruction, the np — > dir data were binned 
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Fig. 10. Comparison of the (d,dp ap ) missing-mass-squared 
distributions at Td = 1.2 GeV when using a polarised hydrogen 
target or filling the cell with nitrogen gas. 
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Fig. 11. Comparison of the (d,pp) missing- mass distributions 
at Td = 1.2 GeV beam energy when using a polarised hydrogen 
target or filling the cell with nitrogen gas. 



in deuteron cm angles. The background subtraction was 
performed separately for each bin and distributions in 
cos(j> built for both spin-up and spin-down modes. The 
weighted sum of the two data sets from different target 
polarisations was taken as the unpolarised mode. Weights 
were determined according to the relative asymmetries 
with respect to the unpolarised state. The ratios of the 
difference to the sum of the data for the two polarised 
modes were than fitted with a linear function in cos</> 
and the value of the product QA y deduced. Taking the 
mean analysing power (A y ) in each 6*^ m bin from the SAID 
pp — > dn + database [31], this gave = 0.61 ± 0.02 and 
Q l = -0.70 ±0.03. 

Using the quasi-free np — > dn° reaction, in an anal- 
ogous way to that for the target polarisation, the vector 
polarisation of the deuteron beam (state 3 in Table [T]), 
was determined to be P, = —0.51 ± 0.05. 




150 
q [MeV/c] 



Fig. 12. Proton analysing powers Ay for the dp — > {pp} a n re- 
action at Td — 1.2 (red squares) and 2.27 GeV (blue triangles) 
for E pp < 3 MeV. The error bars include the uncertainties 
from the target polarisation. Curves correspond to the theo- 
retical predictions. Note that at 2.27 GeV the A? prediction is 
very small and hardly visible on this scale. 



6.4 Measurement of the deuteron-proton 
spin-correlation parameters 



The double-polarised experimental data from the 2009 
beam time offers an excellent opportunity for studies of 
the spin-correlation parameters C x>x and C ViV to deter- 
mine the relative phases of the spin-spin amplitudes [20] . 
Furthermore, it is of interest to see whether the suspected 
£ deficiencies in the SAID amplitudes at highest beam en- 
ergy are reflected also in the spin correlations. 

The background subtraction was carried out in the 
same way as for the dp — » dTT°p sp reaction, with the nitro- 
gen gas data simulating the shape of the background, as 
illustrated in Fig. [TT1 

In order to extract the spin-correlation parameters from 
the dp —> {pp} s n reaction, no tensor polarised deuteron 
beam modes were used (see Table [1]). In this case the ra- 
tio of the polarised N(q,<j>) to unpolarised N°(q) yields 



has the form [31] : 

= 1 + QAV{q) cos + lP z A d y {q) cos + 
+|P,Q[(1 + cos2(p)C y ^ y {q) + (1 - cos2^)C XllB (?)].(6.3) 

Although the experiment was designed for the study 
of spin correlations, by analysing first the polarised target 
yields for an unpolarised deuteron beam, we could obtain 
the dependence of the target analysing power A^ on q 
that is presented in Fig. [TJ at 1.2 and 2.27 GeV. This 
is predicted very well in impulse approximation at the 
lower energy but at 2.27 GeV the corresponding prediction 
can hardly be distinguished from the x-axis on this scale. 
Equation (12.41) then suggests that there must be a serious 
problem also with the SAID determination of the spin- 
orbit amplitude -f(q) at 1.135 GeV. 

In impulse approximation A y vanishes [6], which is 
consistent with our measurements at 1.17 GeV [TO], and 
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Fig. 13. The spin-correlation coefficients C X:X and C y , y for 
the dp — > {pp}„n reaction at Td = 1.2 and 2.27 GeV for 
E pp < 3 MeV. The error bars include the uncertainties from 
the beam and target polarisations. The curves are impulse ap- 
proximation predictions; dashed curves at 2.27 GeV correspond 
to \e(q)\ being reduced by 25%. 



this simplifies the determination of C x>x and Cy iV using 
data with polarised beam and polarised target. After bin- 
ning the normalised counts in intervals in q, the cos 2<j> 
dependence in Eq. (|6.3D allowed us to extract the C x , x 
and C VtV coefficients separately. Note, that the resolutions 
in both E pp and q in the cell-target data are similar to 
those achieved with the cluster target. 

The spin-correlation data for an E pp < 3 MeV cut are 
compared with theoretical predictions in Fig. 1131 The good 
agreement with the experimental points at T c i = 1.2 GeV 
shows that the two relative phases between the spin-spin 
amplitudes are well predicted by the SAID program at this 
energy. It can, however, not come as a complete surprise to 
find that there are serious discrepancies at Td = 2.27 GeV 
but, as shown by the dashed line, these largely disappear 
if the SAID e(q) amplitude is scaled uniformly by a factor 
of 0.75, i.e., by the same factor that brought agreement 
for the A xx and A yy observables! 



7 Conclusions and Outlook 

We have measured the unpolarised differential cross sec- 
tion and the Cartesian tensor analysing powers in the 



dp — ► {pp} s n reaction for small momentum transfers be- 
tween the proton and neutron by using a hydrogen cluster 
target in combination with a tensor polarised deuteron 
beam. The cross section data at 1.2, 1.6, and 1.8 GeV 
are very well described in impulse approximation using 
the current SAID solution for the np amplitudes. These 
amplitudes suffer from much bigger ambiguities at higher 
energies and the corresponding cross section prediction is 
about 15% too high compared to our results at 2.27 GeV, 
though one must bear in mind the 6% uncertainty in the 
overall normalisation of the data at this energy. The sus- 
picion must fall on the SAID solution, which predicts an 
unpolarised np — > pn cross section that may be up to 10% 
too large [2J, of which the spin-dependent contribution is 
also about 10% too large @]. 

The description of A xx and A yy is also very good at 
the three lower energies but much poorer at 2.27 GeV. 
Since one would expect the impulse approximation to be- 
come better as the energy is raised, attention is once again 
focussed on the SAID np amplitudes. The strength of 
|e(0)| 2 in np charge exchange relative to |/3(0)| 2 = |<5(0)| 2 
is determined by the spin-transfer parameters Kll(Q) and 
Knn(0), but there are no measurements of these quanti- 
ties in the relevant angular and energy region. This limits 
severely the SAID predictive power for the deuteron ten- 
sor analysing powers. To fit our data, we have reduced 
the SAID prediction for e(q) uniformly by 25% and this 
reproduces the results much better. Although this might 
be improved further by introducing a g-dependence in this 
factor, the present data do not justify such a refinement. 

By replacing the unpolarised hydrogen cluster-jet tar- 
get by a polarised hydrogen gas cell, it was possible to 
measure the spin-correlation coefficients C x . x and C y . y in 
the dp — ¥ {pp} s n reaction with a vector polarised deuteron 
beam, but only at 1.2 and 2.27 GeV. The behaviour seen 
here is similar to that for the other observables, with a 
good description being achieved at 1.2 GeV whereas at 
2.27 GeV a reduction of the order of 25% seems to be 
required in the e input. 

As a by-product of the polarised cell experiment, we 
were also able to measure the proton analysing power in 
the reaction. As with the other observables, impulse ap- 
proximation reproduces well the small A^ signal at 1.2 GeV 
but fails completely at 2.27 GeV. This suggests that the 
SAID solution for the 7-amplitude is also unreliable at the 
higher energy. 

In summary, the fact that the impulse approximation 
with the current SAID input reproduces well all our data 
below 1 GeV per nucleon gives us confidence that the 
charge- exchange methodology works well. However, the 
discrepancies seen at the higher energy can only be re- 
solved by reducing the strength of the spin-spin ampli- 
tudes, especially in the longitudinal direction, while in- 
creasing the spin-orbit contribution. It is therefore evident 
that the charge exchange on the deuteron contains valu- 
able information on the neutron-proton amplitudes. The 
challenge is to get this used inside the SAID program. 

The experiments reported here were carried out up to 
the maximum deuteron energy available at COSY. To go 
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higher in energy at this facility, an experiment would have 
to be undertaken in inverse kinematics with a polarised 
proton incident on a polarised deuterium gas cell |12) , with 
the two slow protons being detected in the Silicon Track- 
ing Telescopes [35] ■ This will allow the studies reported 
here to be continued up to 2.9 GeV per nucleon [9ll5fj]. 

We are grateful to the accelerator crew for the reliable oper- 
ation of COSY and the deuteron polarimeters. We would like 
to thank I. I. Strakovsky for many useful discussions and for 
providing us with values of the current SAID neutron-proton 
amplitudes. The work was supported by the COSY FFE pro- 
gramme and the Shota Rustaveli National Science Foundation 
Grant 09-1024-4-200. 
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